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INTRODUCTION 


To  ensure  uninterrupted  operation  of  critical  underground  power 
and  commuuica -ion  installations  during  and  after  nuclear  attack 
conditions,  it  is  necessary  that  such  facilities  be  isolated  from 
or  extensively  hardened  against  overpressures  and  ground  motion. 

The  extent  of  hardening  of  such  facilities  will  depend  upon  the 
design  concept  of  the  entire  installation,  expected  overpressures, 
the  extent  of  redundancy  to  assure  facility  survival,  availability 
and  feasibility  of  obtaining  equipments  which  can  withstand  attack 
conditions,  costs,  and  a  number  of  other  operational  considerations. 

In  any  extensive  hardening,  it  is  necessary  to  consider  the  possibility 
of  completely  sealing  off  the  entire  installation  from  dependence 
on  any  surface  facility.  A  major  consideration  associated  with  such 
isolation  is  the  problem  of  dissipation  of  the  power  system's  waste 
heat.  The  requirement  for  a  large  heat  sink  becomes  obvious  for 
any  major  installation. 

The  common  method  of  heat  rejection  during  preattack  operation  of 
such  facilities  would  probably  be  to  the  surface  auove  ground  using 
air  or  water  as  the  medium.  This  is  a  necessary  requirement  for 
facilities  which  are  expected  to  be  operating  most  of  the  time 
during  peacetime.  However,  such  installations  may  be  required 
to  function  in  a  closed  cycle  during  and  after  attack  conditions. 

This  requirement  makes  it  necessary  that  the  waste  heat  be  stored 
underground  for  a  specified  period  of  time. 

Several  alternative  methods  of  heat  storage  for  underground 
facilities  have  been  considered  in  prior  studies  for  conditions 
where  the  heat  storage  systems  were  located  so  deep  in  the  ground 
that  they  were  not  expected  to  be  subjected  to  significant  overpres¬ 
sures  or  ground  motion  during  attack.  For  such  systems,  stored 
chilled  water,  ice-water  mixtures,  and  solutions  of  various  chemicals 
in  water  were  analyzed  as  possible  heat  sinks  [1]  .  Detailed 
analyses  of  these  heat  sinks  and  a  consideration  of  rock  heat  sinks 
have  been  carried  out  by  various  U.  S.  Governmental  agencies,  some 
of  which  are  listed  as  references  [2-13]  in  this  report.  The  power 
source  considered  in  these  analyses  was  a  nuclear  reactor  of  capacity 
in  the  megawatt  range.  A  later  experimental  study  [W]  reports 
the  results  of  a  test  program  carried  out  to  determine  the  feasibility 
of  a  rock  heat  sink.  The  results  showed  severe  structural  failures 
due  to  thermal  stresses  and  moisture  penetration  in  the  rocks. 

For  facilities  which  require  power  in  the  order  of  100  to  500  kw 
and  which  are  located  closer  to  the  surface  so  that  they  are 
subjected  to  the  overpressures  and  ground  moLion  during  attack, 
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the  design  considerations  are  entirely  different  from  those  considered 
in  studies  cited  above.  For  such  applications,  although  storage  of 
ice  or  an  ice-water  mixture  is  possible,  it  may  be  very  costly  if 
the  specified  period  of  post  attack  operation  is  more  than  a  few 
days.  Besides  cost,  the  ice  or  ice-water  mixture  heat  sinks  would 
also  require  regular  or  periodic  surface  support  to  remove  the  heat 
gained  from  surrounding  structures  or  ground  during  the  standby 
period.  Another  possible  heat  sink  could  be  the  presence  of  aquifers 
in  the  vicinity  of  the  installation.  To  use  aquifer  as  a  heat  sink, 
shock-isolated  wells  would  have  to  be  drilled  for  withdrawal  and 
reinjection  of  water.  The  parameters  foi  aquifer  development  and 
use  are  highly  site-dependent,  requiring  extensive  geological 
determination  and  evaluation  at  uach  proposed  site  [15-18]  . 

Another  method  of  underground  storage  of  heat  is  the  use  of 
earth  as  a  heat  sink.  If  the  temperature  of  the  soil  surrounding  an 
underground  installation  could  be  raised,  it  can  act  as  an  excellent 
heat  sink.  To  accomplish  this,  excavation  of  the  ground  would  have 
to  be  carried  out  to  permit  installation  of  pipe  grids  followed  by 
burying  the  grids  in  select  backfill  of  high  thermal  conductivity 
and  specific  heat.  The  spacing  of  pipes  in  the  grid,  pipe  sizes, 
dimensions  of  the  earth  heat  sink  including  allowances  around  it  to 
reduce  the  effects  of  ground  motion  to  an  acceptable  minimum,  the 
amount  of  moisture  in  the  backfill,  and  the  density,  particle  size 
and  type  of  backfill  are  among  factors  which  will  require  careful 
design  and  evaluation. 

This  report  begins  by  summarizing  the  efforts  of  a  literature 
survey  carried  out  to  gather  available  information  on  soils  and 
their  properties  which  determine  the  volume  of  soil  that  will  be 
required  to  absorb  a  predetermined  amount  of  heat  over  a  specified 
period  of  time.  To  ascertain  the  practicality  of  earth  heat  sinks  for 
such  applications,  preliminary  calculations  were  then  carried  out 
to  determine  the  order  of  magnitude  of  dimensions  of  the  required 
heat  sink.  These  calculations  have  shown  that  earth  heat  sinks  may 
be  practical  and  that  they  should  be  considered  along  with  other 
methods  in  the  selection  process  for  a  heat  dissipation  medium  for 
installations  operating  on  a  closed  cycle  basis  over  limited 
durations  in  time. 

GENERAL  CONSIDERATIONS 

The  use  of  earth  as  a  heat  sink  or  source  requires  very  different 
methods  of  analysis  from  those  used  when  the  atmosphere  or  a  large 
body  of  liquid  is  used.  For  air,  water  or  ice  heat  sinks,  the  heat 
transfer  does  not  produce  temperature  gradients  in  the  medium  because 
of  convection;  hence,  analytical  methods  involving  steady  state 
analyses  are  adequate.  When  heat  is  added  to  or  extracted  from  earth, 
however,  heat  diffusion  occurs  so  slowly  that  marked  temperature 
gradients  are  established  which  vary  as  a  function  of  time.  For  such 
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earth  embedded  heat  transfer  surfaces,  the  design  must  be  in  terms 
of  transient  heat  flow,  requiring  transient  heat  conduction  analyses. 

The  importance  of  transient  analysis  is  brought  out  by  the  fact  that 
variations  in  heat  transier  rates  with  time  are  of  the  order  of 
thousands  of  percent.  For  example,  if  a  large  flat  plate  buried  in 
earth  is  raised  suddenly  by  35°  F  above  the  steady  state  earth 
temperature  and  maintained  at  that  level,  the  rate  of  heat  loss 
from  the  flat  plate  after  1,  16  and  3600  hours  will  be  approximately 
20,  5  and  1/3  Btu  per  hour  per  square  foot  [19]  , respectively, 
for  a  soil  of  diffusivity  of  0,025  square  ft/hr.  Thus,  it  is  very 
evident  that  any  statement  as  to  the.  heat  transfer  rate  from  buried 
surfaces  is  meaningful  only  if  the  heat  transfer  rate  is  an  instantane¬ 
ous  one  applicable  after  a  specific  number  of  hours  subsequent  to 
initial  operation. 

Any  underground  power  source  which  is  used  only  under  attack  or 
post  attack  conditions,  when  the  peacetime  surface  power  sources  are 
out  of  operation,  may  be  designed  to  be  used  only  once,  and  that  too 
for  a  limited  duration  in  time  of  perhaps  not  more  than  two  weeks  or 
so.  For  such  a  short  duration,  compared  to  the  use  of  earth  as  a 
heat  source  or  sink  for  year-round  operation  in  a  heat  pump  cycle  for 
home  conditioning,  the  heat  transfer  rates  from  buried  pipes  may  be 
reasonably  high  at  the  end  of  the  designed  operating  period  to  offer 
earth  as  a  useful  alternative  heat  sink  for  underground  power  sources. 
This  may  be  especially  true  if  care  is  taken  in  selecting  a  soil 
of  high  thermal  conductivity  and  specific  heat. 

Beyond  depths  of  just  a  few  feet,  the  soil  temperature  at  any 
location  has  a  steady  year-round  value  which  remains  largely  unaffected 
by  seasonal  climatic  changes  on  the  surface.  The  steady  state  soil 
temperature  at  a  depth  of  30-60  feet  may  vary  anywhere  from  37°F 
to  72°F  in  the  United  States,  depending  upon  the  location  [20]  . 

Below  a  depth  of  about  30  feet,  the  steady  earth  temperature 
further  increases  by  about  1°F  for  each  60  feet  increase  in  depth. 
Therefore,  for  all  practical  purposes,  the  earth  temperature  at  a 
given  location  may  be  assumed  to  have  a  constant  value  before  the 
heat  exchanger  embedded  in  the  ground  has  been  used  either  to  extract 
heat  from  or  to  add  heat  to  it. 

REVIEW  OF  RELATED  WORK 

During  the  1940' s  and  early  1950's,  earth  was  seriously  considered 
as  a  possible  source  of  storage  medium  for  heat  in  the  United  States. 
Investigations  were  carried  out  into  the  feasibility  and  economics  of 
using  earth  in  a  heat  pump  cycle  for  year-round  air  conditioning  of 
homes.  As  heat  pump  cycles  turned  out  to  be  more  expensive  than  gas  or 
oil  heating  systems,  interest  in  them  subsided  causing  a  virtual  halt 
of  investigations  in  this  area.  Most  of  the  published  work  in  heat  pump 
investigations  is,  therefore,  20-30  years  old.  In  the  following  survey, 
the  properties  of  soils  which  determine  the  transient  heat  transfer 
rates  will  first  be  discussed,  followed  by  the  analytical  considerations 
and  test  results  as  presented  in  the  literature. 
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Soil  Properties 


The  soil  properties  which  determine  the  rate  of  heat  transfer 
between  an  embedded  pipe  and  surrounding  earth  are  density,  porosity, 
fluid  permeability,  moisture  contents,  thermal  conductivity,  specific 
heat,  thermal  diffusivity  and  particle  size.  Properties  like  specific 
heat  and  moisture  contents  are  determined  by  the  chemical  and 
mineralogical  composition  of  the  soil,  water  tables  and  porosity. 
Properties  like  thermal  conductivity,  porosity  and  permeability  are 
strongly  dependent  upon  the  form  of  rock,  particle  size  and  packing 
arrangement.  All  of  these  properties  affect  heat  transfer  rates 
through  their  affect  upon  the  thermal  diffusivity  of  the  soil. 

The  constituents  of  soil  vary  from  place  to  place  as  do  particle 
sizes  and  distribution,  density,  porosity  and  the  amount  of  moisture 
present.  It  is,  therefore,  obvious  that  the  thermal  diffusivity  of 
soils  will  vary  widely  from  one  location  to  another.  To  provide 
some  idea  of  these  variations,  the  various  types  of  rocks  subjected 
to  thermal  conductivity  determination  in  one  investigation  may  be 
cited  [ 21-22]  .  Some  of  the  rock  samples  studied  in  this  effort 
were:  granite,  monzonite,  tonalite,  augite,  syenite,  albitite, 
anorthosite,  diabase,  gabbro,  hypersthenite,  bronzitite,  dunite, 
sandstone,  marble,  gneiss,  slate,  limestone,  dolomite,  calcite, 
halite,  and  quartz.  Each  of  these  samples  contained  varying  per¬ 
centages  of  other  minerals,  and  the  crystal  size  varied  between 
0.001  and  3.0  millimeters.  The  thermal  conductivity  of  these  un¬ 
fractured  dry  rocks  was  found  to  vary  between  1.0  and  2.7  Btu/hr-f t-°F. 
The  density  variation  was,  however,  smaller  since  these  materials 
were  in  solid,  dry  form.  In  powdered  form,  with  particles  of  size 
between  10  and  500p,,  the  density  of  inorganic  soils  varies  between 
60  and  120  lbs/'cu  ft.  Wide  variations  in  soil  composition  in  one 
state,  and  significant  variations  with  depth  at  individual  sample 
locations  have  also  been  reported  [23]  .  Some  other  published 
values  of  thermal  conductivity  of  unfractured  rocks  are:  0.7  to  1.7  j 
Btu/hr-f t-°F,  [24]  ,  1.9  Btu/hr-ft°F,  [25]  ,  and  1.52  to  3.46 
Btu/hr-f t-°F,  [26]  A  summary  of  other  published  values  of  the 
thermal  conductivity  of  solid  rocks  ranging  from  0.6  to  3.5  Btu/hr-ft-°F 
is  given  in  reference  [27]  . 

For  the  purpose  of  earth  heat  sink  applications,  however,  the 
important  thermal  conductivity  and  diffusivity  values  are  those  of 
soils  in  powdered  form  containing  varying  percentages  of  various 
minerals  and  water.  If  the  water  tables  are  high,  the  percentage 
of  moisture  will  depend  largely  upon  the  porosity  of  soil.  A  powder 
of  uniform  size  solid  spheres  would  have  a  porosity  between  25.957. 
(rhombohedral  packing  arrangement)  and  47.647.  (cubic  packing  arrange¬ 
ment).  Since  soil  particles  at  any  location  are  likely  to  consist 
of  irregular  grain  shapes  of  sizes  varying  over  a  wide  range,  a 
much  larger  variation  in  porosity  should  be  expected.  Furthermore, 


4 


a  loosely  sift.ed  powder  will  have  a  much  larger  porosity  than  one 
which  has  been  compacted  by  mechanical  means.  Porosity  variations 
for  powder  and  granular  materials  have  been  found  to  vary  from  33-507. 
for  Olivine  Basalt  Powder  to  63-757.  for  Chondrite  powder  [27]  . 

Several  experimental  investigations  to  determine  the  thermal 
conductivity  of  soils  containing  known  amounts  of  moisture  have  been 
reported  in  the  literature  [28-30]  .  Smith  and  Yamauchi  [28]  tested 
five  different  soils,  each  with  several  different  moisture  concentrations. 
Each  soil  was  first  examined  to  determine  the  distribution  of  particle 
sizes.  Of  five  samples,  three  were  classed  as  Sandy  Type  Soils  with 
grain  sizes  of  between  0.08  and  4.0  millimeter  for  707.  of  the  soil.  The 
other  two  soils  were  Clay  and  Loam  with  grain  sizes  of  between  0.005 
and  0.08  millimeter  for  807.  of  the  soil.  The  moisture  contents  as  a 
percentage  of  dry  weight,  densities  and  thermal  conductivities  at 
room  temperature  for  the  two  types  of  samples  are  given  in  Table  1 
to  provide  some  idea  of  the  effects  of  moisture  and  density  upc. a 
the  thermal  conductivity.  Detailed  curves  are  provided  in  the  original 
publication  which  show  strong  dependence  of  thermal  conductivity  of 
a  particular  type  of  soil  upon  both  the  percentage  of  moisture 
present  in  the  soil  and  the  density  of  the  soil.  A  more  detailed 
investigation  of  thermal  conductivities  of  soils  carried  out  for  the 
U.  S.  Army  Corps  of  Engineers  is  available  [  29  ]  .  Properties  of  13 
different  soils,  each  with  several  different  moisture  concentrations, 
are  summarized  in  reference  30. 

An  examination  of  published  literature  suggests  that  the  thermal 
conductivity  of  soil  may  vary  between  around  0.3  and  1.5  Btu/hr-f t-°F, 
depending  upon  the  composition  of  the  soil,  grain  size  and  the  amount 
of  moisture.  If  select  backfills  were  used  for  the  purpose  of 
providing  earth  heat  sinks,  it  is  quite  easy  to  select  soils  with  high 
thermal  conductivities  and  specific  heats.  However,  the  percentage  of 
moisture  will  depend  to  some  extent  upon  the  su/rounding  ground  or 
structure.  If  the  ground  is  relatively  dry,  water  from  the  moist, 
select  backfill  will  diffuse  into  it  with  time  thereby  reducing  the 
effectiveness  of  the  heat  sink.  For  the  purpose  of  preliminary 
calculations  to  determine  the  order  of  magnitude  of  the  heat  sink 
dimensions,  a  conservative  thermal  conductivity  value  of  0.5  Btu/hr-f t-°F 
may  be  assumed. 

The  other  two  properties  of  the  soil  which  determine  heat  transfer 
rates  under  transient  conditions  are  the  density  and  the  specific  heat 
of  the  soil.  The  density  of  inorganic  soils  vary  between  60  and  120 
lbs/cu  ft,  depending  upon  their  composition.  For  the  purpose  of 
preliminary  analysis,  a  soil  density  of  90  lbs/cu  ft  appears  reasonable. 

The  specific  heat  of  the  soil  will  depend  strongly  upon  the 
percentage  of  moisture  present  since  the  specific  heat  of  water  at 
room  temperature  of  1.0  3tu/lb-°F  is  about  four  to  five  times  that  of 
rocks.  For  moist  soils,  values  of  the  specific  heat  considered  in 
the  literature  [l9,  31,  32]  range  between  0.20  and  0.45  Btu/lb-°F.  A 
reasonable  value  for  this  perliminary  analysis  may,  therefore,  be 
assumed  to  be  0.3  Btu/lb-°F. 


5 


With  assumed  values  of  the  thermal  conductivity  of  0.5  Btu/hr-f t-°F, 
density  of  90  lbs/cu  ft  and  specific  heat  of  0.3  Btu/lb-°F,  the  value 
of  soil  thermal  diffusivity  for  sample  calculations  is  0.0185  square 
ft/hr. 


Analytical  Considerations 


Consider  an  infinitely  long  permanent  line  source  or  sink  for 
heat  in  an  infinite  soil  medium  at  an  initial  uniform  temperature, 
i.e.,  constant  heat  flux  at  the  pipe  wall.  Subsequent  temperatures 
at  any  point  in  the  medium  may  be  represented 

I  =  Io  +  25k/T-dP 

A 


(i) 


where  Tq  =  Initial  uniform  temperature  of  soil,  °F. 

T  =  Temperature  of  soil  at  time  t  and  distance  r,  °F. 

Q'  =  Heat  Transfer  rate  across  the  pipe,  Btu/hr  per  foot  length, 
k  =  Thermal  conductivity  of  soil,  Btu/hr-f t-°F . 

X  2 

r  =  Distance  from  the  center  of  pipe,  ft. 
a  =  Thermal  diffusivity  of  soil,  k/pCp  ,  sq  ft/hr. 
p  =  Density  of  soil,  lb/cu  ft. 

=  Specific  heat  of  soil,  Btu/lb-°F 
t  =  Time  from  start  of  operation,  hr. 

P  =  Variable  of  integration. 


Ingersoll  and  Plass  [32]  offered  Equation  (1)  in  the  form 

AT  =  I(X) 


where 


I(X)  = 


l 


00  -p 


dp 
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and  provided  values  of  the  integral  I(X)  for  values  of  X  ranging 
between  0.0001  and  2.20  in  tabular  form. 

Although  Equation  (2)  is  exact  only  for  a  true  line  source,  it 
may  be  applied  with  negligible  error,  after  a  few  hours  of  operation, 
to  small  diameter  pipes  in  actual  use  in  heat  pump  installations. 

In  fact,  for  Fourier  numbers  (N^)  of  <.  0.05,  the  error  [35]  in  using 

Equation  (2)  is  less  than  2 7„.  For  values  of  the  Fourier  number 
2 

(r  /at)  larger  than  0.05,  a  more  general  equation  which  is  applicable 
for  all  values  of  the  Fourier  number  may  be  used.  This  general 
equation  for  a  constant  heat  transfer  rate  may  be  written  in  the 
form  [36]  , 


T  -  T 

o 


JpCpP^CP)  -  Jt  (P)Y0(pP) 
J*(P)  +  y12(P) 


d£  (3) 

„2 


Where  Z  =  =  & 

F  r 

P  =  r/rL 

r^  *  Radius  of  the  pipe 


Jq,  =  Bessel  functions  of  the  first  kind  of  order  zero  and 
one,  repsectively . 


Yq,  Y1  =  Bessel  functions  of  the  second  kind  of  order  zero  and 
one,  respectively. 


For  earth  heat  sinks  for  underground  applications,  the  pipe 
length  is  going  to  be  several  hundred  times  the  pipe  diameter. 
Furthermore,  as  shown  later,  the  values  of  Z  are  going  to  be  much 
larger  than  20.  Therefore,  Equation  (2)  may  be  adequate  for 
applications  using  single  pipes  buried  in  the  ground. 

It  is  very  unlikely,  however,  that  a  single  pipe  will  be 
adequate  to  handle  heat  loads  of  more  than  just  a  few  thousand  Btu's 
per  hour,  unless  the  pipe  is  made  several  miles  long.  From  a 
realistic  and  cost  effective  viewpoint,  however,  it  will  be  necessary 
to  use  several  "hairpin  loop'  pipes  in  parallel.  Ingersoll  and  Plass 
[  32]  considered  two  parallel  pipes  of  radius  r^,  r ^  feet  apart,  in  a 
hairpin  loop,  and  proposed 
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Where 


T.  -  T  = 
1  o 


*  IW  [I(xi 


)  +  I(X,) 


=  r1/2  Vat 
=  r2/2  Vat 


T^  =  Temperature  at  pipe  wall  of  radius  r^. 

For  an  actual  ear  .  heat  exchanger,  practical  considerations  would 
dictate  the  use  of  long  parallel  pipes  in  a  3-dimensional  array.  For 
such  an  application,  there  is  no  reason  why  Equatior  (4)  cannot  be 
mofidied  to 


T.  -  T 
1  o 


_  s  Kxn) 

27Tk  n=l  n 


Where  X  =  r  /2  V"^t 
n  n  u 

and  r2,  r^,  r^,  etc.,  are  distances  of  surrounding  pipes  in  the  array 
from  the  center  of  the  pipe  of  radius  r^  under  consideration. 

A  somewhat  different  equation  will  result  if  instead  of  maintaining 
a  constant  heat  flux  on  the  surface  of  the  pipe,  a  constant  pipe  wall 
temperature  was  maintained.  For  a  buried  pipe  whose  temperature  is 
raised  above  or  lowered  below  the  uniform  earth  temperature  and 
maintained  at  that  level,  the  heat  transfer  rate  will  decrease 
asymptotically  with  time.  Equation  for  this  constant  pipe  temperature 
case  may  be  written  in  the  form  [37]  : 


=  8KAT  / 

'  »  J  ~j 


02(P>  +  y/us) 


This  equation  has  been  evaluated  by  Ingersoll,  et.al.  [35]  in  the  form 
of  tabulated  values  of  F(Z),  where 


F(Z)  = 


oo 

; 


02<w  +  yo2(P) 
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for  values  of  Z  ranging  between  0.01  and  25,000. 

A  constant  heat  transfer  rate  from  the  surface  of  a  very  long 
cylinder  (thereby  eliminating  end  effects)  embedded  vertically  in  soil 
was  considered  by  Guernsey,  Betz  and  Skau  [3l]  .  Earth  of  infinite 
extent  at  a  uniform  temperature  was  assumed  and  the  effects  of  heat 
interchange  with  atmosphere  was  neglected.  With  these  assumptions, 
they  obtained  an  approximate  solution, 

T  -  To  •  5^75k  [lo810  Z  +  °-78/Z  +  °'351]  ■  (8) 


with  an  error  of  less  than  1%  for  Z  >6.  Using  this  equation,  they 

computed  the  lengths  of  various  diameter  pipes  required  to  extract 

heat  at  a  fixed  rate  over  a  20-year  period  for  a  maximum  predetermined 

value  of  (T  -  T  ) .  They  found  that  increasing  the  pipe  size  from 
o 

3/4-inch  to  3.82-inch,  a  5.1- fold  increase,  decreases  the  required 
length  by  only  20%.  It  therefore  appears  that  smaller  diameter  pipes 
are  more  effective.  This  conclusion  is  reasonable  since  a  longer, 
smaller  diameter  pipe  is  exposed  to  a  larger  volume  of  earth  with 
which  heat  can  be  exchanged. 

All  the  above  equations  ignore  the  effects  of  convection  in  the 
air  or  in  the  moisture  present  in  the  soil  because  it  does  not 
appear  that  convection  currents  in  soils  can  be  significant.  Published 
geological  studies  indicate  that  convection  in  porous  earth  strata  is 
extraordinarily  slow  M  • 

A  somewhat  different  application  has  been  analyzed  by  Kingston  [39] 
which  may  be  important  in  consideration  of  earth  heat  sinks  for  under¬ 
ground  applications.  He  analyzed  the  cooling  of  mass  concrete  by 
water  flowing  in  embedded  pipes.  This  problem  occurs  in  cooling  of 
dams.  He  assumed  that  the  water  flowing  through  a  pipe  cools  the 
concrete  in  a  prism,  the  cross-section  of  which  is  a  hexagon.  The 
hexagon  was  approximated  by  a  circle  and  it  was  assumed  that  the 
outer  boundary  of  the  infinitely  long  cylinder  was  thermally  insulated. 
The  entire  concrete  was  assumed  to  be  initially  at  a  uniform  temperature 
and  subjected  to  a  zero  temperature  at  the  surface  of  a  cylindrical 
hole  (equivalent  to  the  pipe  wall)  concentric  with  the  axis.  The 
differential  equation  to  be  solved  in  this  case  was 


2SL.  = 

dT 


a 


ii  + 
2 
3* 


l  ££ 

r  ^r 


(9) 


with  the  boundary  conditions 


T 


dT 

dr 

Where  r. 

1 


r 

o 


0  at  r  =  r^ 

0  at  r  =  r 

o 

Inner  radius  of  cylinder 
Outer  radius  of  cylinder 


The  solution  of  Equation  (9)  was  quite  involved  and  required  determination 
by  numerical  methods. 

Experimental  Considerations 

Smith  [40]  carried  out  an  experimental  investigation  to  determine 
the  heat  transfer  rates  as  a  function  of  time  for  a  single  isolated 
pipe,  for  two  parallel  pipes  with  varying  j pacings  and  for  three 
parallel  pipes.  He  compared  his  single  pipe  data  with  Equations  (2) 
and  (6)  by  plotting  Q'/k/^T  against  Z  for  both  equations  and  the  data, 
and  found  good  agreement.  The  difference  between  the  two  equations  was 
found  to  be  small,  and  this  difference  was  found  to  decrease  with  an 
increase  in  the  value  of  Z.  The  heat  coefficients  and  the  temperature 
differences  between  pipe  wall  and  initial  earth  temperature  in  the 
test  using  three  pipes  in  parallel  were  found  to  agree  very  well  with 
values  predicted  by  Equation  (5).  The  effects  of  the  presence  of 
other  pipes  upon  the  heat  transfer  rates  from  a  pipe  were  presented  in 
graphical  form.  The  ratios  of  heat  transfer  rates  from  a  pipe  in  arrays  of 
of  three  and  infinite  number  of  pipes  to  its  value  for  a  single  pipe 
were  also  plotted  against  Z  for  various  spacings  between  pipes  in  the  arrays. 

Hadley  [4l]  correlated  all  available  experimental  data  from  heat 
pump  ground  coil  installations  from  eight  different  sites  by  plotting 
K/^T/Q1  against  Z  and  found  that 

^  =  0.311  log1QZ  -  0.3696  .  (10) 

When  he  compared  Equation  (10)  with  Equation  (2),  he  found  that 
Equation  (2)  tends  to  predict  conservative  values  for  design  purposes. 

For  example,  the  experimentally  obtained  value  of  Q'/k^T  from  Equation  (10) 
was  about  2.25  times  the  value  predicted  by  Equation  (2)  at  Z  =  100. 

At  Z  =  1000,  the  experimental  value  of  Q'/k/^  T  was  approximately  307. 
hig’er  than  the  predicted  value.  The  difference  approached  zero  as  Z 
approached  20,000.  From  this  comparison,  it  may  be  readily  seen  that 
Equation  (2)  should  provide  somewhat  conservative  designs. 
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THEORETICAL  CONSIDERATIONS 


The  use  of  an  earth  heat  sink  to  absorb  large  quantities  of  heat 
over  a  relatively  short  period  of  time  requires  a  large  surface 
area  of  the  pipe,  but  a  small  soil  volume  because  heal  transmitted 
to  soil  at  the  pipe  wall  can  diffuse  only  a  short  distance  in  the  soil 
over  the  anticipated  2-3  weeks  of  operation.  To  minimize  excavation 
and  backfill  costs,  therefore,  it  seems  reasonable  to  install  long 
pipes  in  such  a  way  as  to  enclose  the  required  amount  of  soil  in  a 
parallelopiped  arrangement.  Either  of  the  two  pipe  arrangements 
whose  cross-sections  are  shown  in  Figures  1  and  2  is  adequate  for 
this  application.  For  the  purpose  of  thermal  analysis,  both  of  these 
arrangements  could  be  approximated  by  a  bank  of  hollow  soil  cylinders 
of  inner  and  outer  radii  of  r^  and  r  .  The  outer  boundary  of  the 

cylinder  at  r  =  r^  may  then  be  assumed  to  be  perfectly  insulated  so 

that  no  heat  transfer  can  take  place  across  its  boundary.  Although 
this  assumption  is  incorrect  for  cylinders  which  form  the  boundaries 
of  the  matrix,  their  effect  upon  the  total  heat  transfer  capability 
of  the  matrix  may  be  assumed  to  be  negligible  if  the  matrix  has  a 
large  number  of  pipes  and  if  the  dimension  r^  is  such  that  the  heat 

transfer  rate  from  pipes  of  radius  r ^  after  the  specified  time 

differs  little  from  that  of  the  pipe  located  alone  in  an  infinite  earth 
medium.  Furthermore,  if  the  pipes  are  long  enough,  the  end 
effects  may  be  assumed  to  be  negligible.  The  hexagonal  prism 
arrangement  of  Figure  2  provides  a  better  approximation  becuase  it 
encloses  only  9.31%  of  earth  volume  that  will  not  be  accounted  for 
in  this  analysis  compared  with  the  unaccountable  21.46 %  in  the 
Figure  1  arrangement.  With  these  assumptions,  the  appropriate 
unsteady  state  heat  transfer  equation  that  has  to  be  solved  is  the 
same  as  Equation  (9),  but  with  different  initial  and  boundary 
conditions.  To  simplify  analysis,  the  zero  point  of  the  scale  can 
be  shifted  to  the  steady  state  initial  earth  temperature.  The 
temperatures  thus  obtained  from  the  solution  of  Equation  (9)  will 
represent  the  rise  of  temperature  above  the  steady  state  earth 
temperature.  The  boundary  conditions  are,  therefore, 

T  =  T,  at  r  =  r. 

1  i 


at  r 


r 

o 


The  solution  of  this  equation  involves  Bessel  functions  of  the 
first  and  second  kind  and  of  the  first  and  second  order.  The  solution 
is,  therefore,  quite  involved,  time  consuming,  and  requires  the 
extensive  use  of  a  high  speed  computer.  This  will  be  carried  out  for 


wide  ranges  of  ar>  Z  and  t  during  the  second  phase  of  this  investigation. 
For  the  purpose  of  this  report,  however,  preliminary  estimates  of 
the  sizes  of  earth  heat  sinks  were  carried  out  as  described  below  using 
Equation  (5). 


PRELIMINARY  CALCULATIONS 
Design  Criterion 

Inexpensive  power  sources  in  the  100  kw  to  500  kw  range  are  non¬ 
nuclear  in  nature.  Suitable  non-nuclear  power  sources  that  might  be 
considered  for  underground  applications  include  Diesel,  Gas  Turbine 
and  Fuel  Cells.  The  combustion  products  or  exhausts  of  Diesel 
power  sources,  for  example,  are  at  temperatures  approaching  1000°F. 

Since  the  efficiency  of  these  power  sources  does  not  depend  upon  the 
sink  temperature,  it  should  be  possible  to  operate  the  heat  sink  at 
fairly  high  temperatures.  However,  because  of  the  recycle  nature 
of  closed  cycle  operation,  the  operating  temperature  limit  of  the 
earth  heat  sink  may  be  determined  by  the  temperature  of  the  processed 
exhaust  which  has  to  be  recycled  to  the  power  source.  For  the 
purpose  of  this  discussion,  only  a  Diesel  system  will  be  considered. 

In  a  closedcycle  diesel  operation,  the  inlet  temperature  of  fuel, 
oxidizer  and  diluent  (or  heat  absorber)  has  to  be  about  180°F  if  the 
engine  is  operating  in  a  psychrocycle  mode,  to  ensure  the  presence  of 
sufficient  moisture  to  ptevent  the  cylinder  temperature  from  rising 
beyond  design  limits.  If  the  engine  is  to  operate  on  a  kreislauf 
cycle,  however,  the  inlet  temperature  has  to  be  somewhat  higher  than 
this  value  to  ensure  combustion.  In  either  case,  the  engine  cooling 
jacket  or  radiator  temperature  is  limited  to  about  212°F  when  using 
ebullient  cooling  in  commercially  available  units.  The  earth  heat 
sink  can  therefore  be  subjected  to  a  step  temperature  rise  at  pipe 
walls  to  at  least  212°F,  and  possibly  higher.  One  method  of  using 
the  earth  heat  sink  to  absorb  heat  at  212°F  will  be  to  cool  the  combustion 
products  (CO^,  C>2  and  ^0)  from  around  1000°F  to  around  212°F  in  a 

series  of  water  cooled  heat  exchangers  prior  to  the  absorption  of  carbon 
dioxide.  The  steam  thus  formed  in  the  radiator  and  the  heat  exchangers 
can  be  used  to  carry  heat  to  the  earth  heat  sink  and  returned  to  the 
radiator  and  the  heat  exchangers  as  liquid  water  at  about  212°F.  The 
pipe  wall  temperature  may,  therefore,  be  assumed  to  be  subjected  to 
212°F  over  the  entire  duration  of  operation. 

The  actual  situation  in  such  a  heat  exchanger  will  be  more  compli¬ 
cated  since,  although  the  temperature  of  steam  entering  and  that  of 
water  leaving  the  heat  sink  has  been  assumed  to  be  at  a  constant 
value  of  212°F,  the  heat  output  rate  of  power  unit  that  has  to  be 
absorbed  by  the  heat  sink  is  also  a  constant.  Therefore,  the  assumed 
initial  and  boundary  conditions  will  cause  a  progressively  larger 
length  of  the  pipe  to  become  active  in  the  heat  exchange  process,  as 
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Che  heat  transfer  rate  from  the  pipe  section  exposed  to  steam  drops 
with  time.  The  heat  transfer  rate  from  any  given  segment  of  the  pipe 
will,  therefore,  vary  with  time  as  will  its  wall  temperature  as  cool 
water  is  replaced  by  warmer  water  which  is  finally  replaced  by  steam 
at  212°F.  The  upstream  sections  of  the  pipe  will  be  the  first  to 
be  exposed  to  212°F  steam  whereas  the  downstream  end  will  be  exposed 
last  at  a  time  when  although  the  entire  heat  exchanger  upstream  is 
being  maintained  at  212°F,  the  total  heat  transfer  rate  is  less  than 
the  total  waste  heat  output  rate  of  the  power  unit.  Similarly,  the 
condensa'.t  temperature  will  probably  start  from  slightly  above  the 
initial  earth  temperature,  rising  to  212°F  at  the  end  of  the  period. 

Obviously,  the  constant  heat  flux  assumption  of  Equation  (5)  is 
not  applicable  for  any  small  section  of  the  pipe,  although  it  is 
applicable  to  the  earth  heat  sink  as  a  whole.  The  constant  pipe  wall 
temperature  case  may  be  somewhat  more  appropriate  here,  although  the 
duration  of  exposure  to  steam  will  vary  from  100%  of  engine  operating 
period  at  entrance  to  the  heat  sink  to  only  a  fraction,  as  yet  undefined, 
at  the  heat  sink  outlet.  For  several  reasons,  however,  Equation  (5) 
was  selected  for  the  preliminary  calculations  shown  here:  (a)  the 
constant  temperature  case,  Equation  (6),  cannot  be  utilized  for  a 
3-dimensional  matrix;  (b)  as  shown  by  Smith  [38]  ,  there  is  only  a 
small  difference  in  heat  transfer  rates  between  the  constant  wall 
temperature  and  the  constant  heat  transfer  rate  cases  for  values  of 
Z  of  500  or  more;  (c)  Equation  (5)  can  be  used  for  a  bank  of  pipes, 
as  confirmed  bv  experimental  investigation  of  Smith  [38]  ,  and 
(d)  the  actual  heat  sink  will  experience  conditions  somewhere  between 
the  constant  temperature  and  the  constant  heat  transfer  rate  cases, 
for  which  equations  are  not  available.  To  permit  application  of 
Equation  (5),  it  was  assumed  that  the  entire  pipe  matrix  is  raised 
to  212°F  at  the  start  of  the  power  unit.  The  heat  transfer  rate 
calculated  at  the  end  of  the  specified  period  was  then  assumed  to 
be  applicable  to  the  whole  pipe  matrix.  This  assumption  will 
obviously  provide  conservative  values  of  the  heat  transfer  rates  since, 
after  the  specified  period,  the  pipe  sections  downstream  will  have  a 
higher  heat  transfer  rate  because  they  have  not  been  exposed  to  212°F 
for  the  entire  test  duration. 

Heat  Transfer  Rate  at  End  of  Specified  Period 

From  Equation  (5),  factors  which  determine  the  heat  transfer 
rate  at  the  end  of  the  specified  period  are:  thermal  diffusivity  and 
the  thermal  conductivity  of  soil,  the  radius  and  spacing  of  pipes,  the 
pipe  wall  to  initial  earth  temperature  difference  and  the  specified 
period  of  time. 

Conservative  values  of  thermal  conductivity  and  thermal 
diffusivity  of  soil  were  earlier  determined  to  be  0.5  Btu/hr-ft-°F  and 
0.0185  sq  ft/hr,  respectively.  A  futher  assumption  may  be  made  that 
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tne  soil  is  at  a  steady  state  temperature  of  60°F  at  the  start  of  its 
use  as  a  heat  sink.  This  assumption  may  again  be  somewhat  conservative 
in  view  of  the  soil  temperature  variation  in  the  United  States  of 
between  37°  and  72°F.  The  value  of  i  T  is  therefore  assumed  to  be 
equal  to  152°F. 

Three  time  periods  will  be  assumed  for  these  calculations:  one 
week,  two  weeks  and  three  weeks.  It  is  doubtful  that  any  emergency 
underground  system  will  be  designed  to  operate  for  much  less  than  one 
week  or  much  more  than  three  weeks.  For  the  purpose  of  preliminary 
calculations,  therefore,  we  have 

k  =  1.5  Btu/hr-f t-°F 

a  =  0.0185  sq  ft/hr, 

AT  =  152°F, 

and  t  =  168,  336  and  504  hours. 

The  last  item  we  need  in  calculating  heat  transfer  rates  from  the 
embedded  pipes  is  the  size  and  spacing  of  pipes.  The  pipe  diameter 
should  not  be  too  small  since  it  may  cause  a  large  pressure  gradient 
in  steam  flowing  through  it,  or  too  large  since  there  is  little  heat 
transfer  improvement  [31]  as  compared  to  the  costs.  Furthermore,  to 
reduce  the  pressure  drop,  several  small  diameter  pipes  may  be  connected 
in  parallel  through  headers.  Based  upon  these  considerations,  a 
pipe  of  outer  diameter  1.5-inch  was  selected.  For  a  pipe  of  radius 
3/4-inch,  the  values  of  Z  for  the  three  selected  time  periods  were 
796,  1591  and  2387,  respectively.  For  these  values  of  Z,  pipe  spacings 
were  determined,  by  trial  and  error,  which  caused  the  heat  transfer 
rates  to  be  approximately  half  those  from  a  single  pipe  located  in 
an  infinite  earth  medium  during  the  same  period.  Acceptance  of  some 
deterioration  in  the  heat  transfer  rate  was  considered  necessary 
to  reduce  the  volume  of  the  earth  heat  sink,  and  the  507>  deterioration 
level  was  chosen  arbitrarily.  Under  the  same  assumptions  and  design 
conditions,  the  pipe  spacing  will  increase  with  an  increase  in  the 
value  of  Z.  In  the  hexagonal  prism  matrix  of  Figure  2,  distances  between 
the  center  of  a  pipe  and  the  centers  of  six  surrounding  pipes,  2rQ, 

were  estimated  to  be  25-inches,  30-inches,  and  40-inches  for  Z  =  796, 

1591  and  2387,  respectively.  Values  of  I(^n)  were  then  obtained  for 

the  different  pipe  locations  in  the  matrix  followed  by  determination  of 
the  average  heat  transfer  rate  as  follows. 

As  a  first  approximation,  24  rows  of  24  pipes  each  were  assumed 
in  the  matrix,  a  section  of  which  is  shown  in  Figure  2.  A  pipe  in 
any  location  was  then  considered  to  be  influenced  by  pipes  no 
further  than  a  distance  of  4r^  away  from  its  center;  i.e.,  only 

including  up  to  two  hexagons  formed  by  pipes  surrounding  it.  This 

assumption  simplified  the  calculations  significantly  whereas  its 

effect  upon  the  heat  transfer  rates  was  found  to  be  small,  as  shown  later. 
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The  number  of  pipes  at  distances  to  4rp  which  affect  a  pipe  location 

are  shown  in  Table  2  for  the  15  different  locations  which  exist  in  the 
24  x  24  matrix. 

For  heat  sources  at  distances  r.,  2r  ,  2  V3~ r  and  4r  .the  values 

o’  o  o’ 

of  I(Xn)  were  calculated  by  linear  interpolation  of  tables  by  Ingersoll  and 

Plass  [32]  for  the  three  time  periods,  1-week,  2-week  and  3-week, 
each  with  its  spacing  of  25,  30  or  40  inches,  respectively,  as  shown  in 
Table  3. 

From  these  values  of  I(X  )  the  value  of  the  heat  transfer  rate 

v  n ' 

for  each  location  was  calculated  from  Equation  (5),  as  shown  in 
Table  4.  It  may  be  noted  from  this  table  that  the  weighted  average 
heat  transfer  rate  for  all  pipes  in  the  24  x  24  matrix  is  69.8  Btu/hr-ft 
which  is  approximately  4%  above  the  value  of  67.1  Btu/hr-ft  for  a  pipe 
at  location  9,  where  it  is  completely  surrounded  by  pipes  on  all  sides 
at  distances  to  4rQ,  for  the  matrix  designed  for  1-week  use.  This 

improvement  due  to  including  the  effect  of  peripheral  pipes  through 
weighted  averaging  is  rather  small  considering  the  extent  of  uncertainty 
in  the  thermal  diffusivity  of  the  soil.  Similar  improvements  of  about 
4%  also  may  be  noted  in  the  calculations  fo  the  2-week  and  3-week  matrix. 

The  effect  of  surrounding  pipes  in  the  matrix  on  the  heat 
transfer  rate  of  a  pipe  may  now  be  evaluated.  Since  about  70%  of  the 
pipes  are  surrounded  by  pipes  on  all  sides  at  distances  to  4rQ, 

similar  to  location  9,  the  effects  of  the  presence  of  pipes  within 

distances  of  2r  and  between  2r  and  4r  on  a  pipe  will  be  considered 
o  00 

separately  for  this  location  for  comparison  purposes.  The  value  of 
I(Xn)  for  a  pipe  located  in  an  infinite  earth  medium  in  the  1-week 

matrix  is  3.778,  Table  3.  If  the  six  surrounding  pipes  at  a  distance 
of  2r  are  considered,  their  effect  is  to  add  6  x  0.397  or  2.382  to  the 

I(Xn)  value  of  3.778.  Similarly,  the  12  additional  pipes  at  distances 

of  over  2rQ,  but  not  greater  than  4rQ  from  location  9  add  0.9552  to 

the  value  of  I(xn).  From  Equation  (5),  therefore,  the  heat  transfer 

rate  for  the  pipe  located  in  an  infinite  earth  medium  is  126.4  Btu/hr-ft 
whereas  the  presence  of  pipes  at  a  distance  of  2rQ  reduces  it  to 

77.5  Btu/hr-ft,  a  reduction  of  39%.  Pipes  located  at  distances  over 

2r  ,  but  up  to  4r  reduce  this  value  further  to  67.1  Btu/hr-ft,  an 
0  0 

additional  reduction  of  about  13% . 

To  determine  the  extent  of  errors  introduced  in  these  calculations 
because  of  not  including  the  effects  of  pipes  at  distances  greater  than 
4rQ,  a  sample  calculation  was  carried  out  for  a  pipe  in  location  16, 

for  which  the  effect  of  an  additional  18  pipe  within  a  radius  of  6r 

0 
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could  be  included,  and  its  heat  transfer  rate  compared  with  that  for  a 
pipe  in  location  9  calculated  earlier.  For  the  1-week  duration  matrix, 
the  heat  tiansfer  rate  for  location  16  was  found  to  be  65.3  Btu/hr-ft, 
a  decrease  of  2.7%  from  the  value  of  67.1  Btu/hr-ft  for  location  9. 

It  may,  therefore,  be  seen  that  placement  of  pipes  in  a  hexagonal 
arrangement  at  distances  of  2rQ,  4rQ  and  6rQ,  as  shown  in  Figure  2, 

causes  successive  reductions  in  the  heat  transfer  rate  of  a  pipe  by 

397o,  13%  and  2.770,  respectively.  The  presence  of  additonal  pipes 

at  distances  over  6r  will  have  an  effect  even  smaller  than  the  2.7% 

o 

caused  by  18  pipes  at  distances  over  4rQ  but  up  to  and  including  6rQ. 

It  may,  therefore,  be  reasonable  to  assume  that  if  the  effect  of  all 
pipes  upon  all  other  pipes  in  the  matrix  were  included,  the  additional 
deterioration  in  the  weighted  average  heat  transfer  rate,  Table  4, 
due  to  pipes  at  distances  greater  than  4r^  may  not  be  greater  than  the 

47»  improvement  in  it  caused  by  peripheral  pipe  locations.  The  average 
heat  transfer  for  the  matrix  may,  therefore,  be  assumed  to  be 
approximately  that  of  the  pipe  in  location  9  when  the  effect  of  pipes 
at  distances  up  to  4rQ  is  included.  Therefore,  in  determining  the  length 

of  pipes  required  for  the  heat  sink,  values  of  Q'  for  the  1-week,  2-week 
and  3-week  durations  will  be  assumed  to  be  67.1,  53.6  and  57.5  Btu/hr-ft, 
respectively.  The  actual  average  heat  transfer  rates  will  be  somewhat 
larger  because  of  the  end  effects  and  because  the  entire  pipe  matrix  is 
not  subjected  to  the  maximum  212°F  over  the  total  time  periods  of  one, 
two  and  three  weeks. 

Dimensions  of  the  Earth  Heat  Sink 

With  a  knowledge  of  the  average  heat  transfer  rate  from  pipes  in 
the  matrix  at  the  end  of  specified  period  of  time,  the  length  of  the 
heat  sink  can  be  determined  from  the  value  of  the  total  heat  that  has  to 
be  dissipated  per  hour  of  engine  operation.  The  excess  heat  generated 
during  a  closed-cycle  operation  of  a  Diesel  engine  is  approximately 
11,500  Btu/hr  per  kilowatt  of  output  in  the  100  to  500  kw  range  of 
engines.  For  the  preliminary  calculations,  an  engine  in  the  midrange 
of  300  kw  may  be  considered  for  durations  of  one,  two  and  three  weeks. 

The  heat  to  be  dissipated  is,  therefore,  at  the  approximate  rate  of 
3.45  x  106  Btu/hr.  From  values  of  Q'  of  67.1,  53.6  and  57.5  Btu/hr-ft 
obtained  earlier  for  an  assumed  24  x  24  pipe  matrix  for  one,  two  and 
three  weeks  of  operation,  the  lengths  of  pipes  required  are  51,400, 

64,300  and  60,000  feet,  respectively.  The  overall  dimensions  of  suit¬ 
able  earth  heat  sinks  are  given  in  Table  5. 

Although  the  overall  volume  of  earth  required  to  absorb  the  waste 
heat  is  large,  the  costs  associated  with  an  earth  heat  sink  are  limited 
to  cost  of  excavation  and  select  backfill  which  is  estimated  at  around 
50<?  per  cubic  foot  [42]  ,  and  the  cost  of  pipes  in  the  matrix.  Furthermore, 
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the  pipe  lengths  and  earth  volumes  calculated  above  were  not  optimized 
for  the  lowest  cost.  An  optimization  involving  tradeoff  between  the 
cost  of  pipes  and  the  cost  excavation  and  backfill  may  reduce  the  total 
cost  of  the  earth  heat  sink. 

For  example,  the  values  of  thermal  conductivity,  specific  heat 
and  soil  density  assumed  in  above  calculations  were  conservative.  If 
care  is  exercised  in  selecting  a  soil  of  high  thermal  conductivity 
(e.g.,  1.2  Btu/hr-f t-°F) ,  high  specific  heat  (e.g.,  0.45  Btu/lb-°F) 
and  a  density  of  100  lb/cu  ft,  then  the  value  of  Q'  for  location  9 
in  the  matrix  using  the  same  size  of  pipes  and  same  spacing  between 
pipes  will  be  approximately  127.6,  104.0  and  111.6  Gtu/hr-ft  for 
operations  of  one,  two  and,  three  weeks,  respectively,  as  shown  in 
Table  6.  The  required  length  of  pipes  may  then  be  reduced  to  27,000, 
33,200,  and  30,900  feet  for  operation  over  one,  two,  and  three  weeks, 
respectively,  and  the  overall  dimensions  of  the  heat  sink  reduced  to 
those  shown  in  Table  7. 

As  may  be  seen,  the  overall  volume  of  the  carefully  selected  soil 
is  approximately  half  that  of  the  soil  with  conservative  thermal 
conductivity  and  specific  heat  values  under  identical  operating 
conditions,  pipe  size  and  pipe  spacing.  If  the  spacing  between  pipes 
were  increased  because  in  a  highly  conducting  soil  the  heat  can 
diffuse  further  during  the  same  time  period,  a  saving  in  the  overall 
lengths  of  the  pipes  of  approximately  5000,  8700,  and  6400  feet  for 
durations  of  one,  two,  and  three  weeks,  respectively,  may  be  obtained 
at  the  expense  of  soil  volume,  as  shown  in  Tables  8  and  9. 

Comparing  the  heat  sinks  with  alternative  spacings,  Tables  7  and  9, 
the  increase  in  the  required  soil  volume  is  20,000,  74, 000, and  90,000 
cubic  feet,  which  is  approximately  equivalent  to  4.0,  8.5,  and  14.0 
cubic  feet  of  soil  per  foot  of  pipe  length  saved  for  the  1-week, 

2-week,  and  3-week  durations,  respectively.  At  an  estimated  cost  of 
50<:  per  cubic  foot  of  pipe  for  excavation  and  backfill,  the  alternative 
heat  sink  of  Table  9  for  the  1-week  service  will  be  more  economical 
if  the  cost  of  the  heat  exchanger  in  the  matrix  is  equivalent  to  less 
than  $2.0  per  foot  of  pipe  used  in  its  construction. 

DISCUSSION 

The  cost  of  an  earth  heat  sink  may  now  be  estimated  to  provide 
a  basis  for  comparison  with  other  heat  dissipation  methods.  A  heat 
sink  matrix  made  from  suitable  stainless  steel  pipes  of  1.5-inch 
diameter  which  can  withstand  external  pressures  of  about  1000  psi 
may  be  assumed  to  cost  approximately  $2.0  per  foot  length  of  pipe. 
Similarly,  a  cost  of  $3.00  per  foot  of  pipe  may  be  assumed  for  service 
at  overpressures  of  about  1500  psi.  The  cost  of  excavation  and  select 
backfill  may  be  assumed  to  be  approximately  50c  per  cubic  foot.  For 
the  purpose  of  this  estimate,  the  costs  associated  with  hardening  the 
heat  sink  against  ground  motion  will  not  be  included.  Using  values 
developed  in  Table  7,  the  cost  of  the  earth  heat  sink  for  1-week  service 


is  $109,000;  the  cost  of  pipes  and  the  cost  of  excavation  and  select 
backfill  being  approximately  equal.  The  total  heat  dissipated  by 
this  heat  sink  over  a  period  of  one  week  will  be  5.79  x  10^  Bta,  or 
approximately  5300  Btu  per  dollar  cost  of  the  earth  heat  sink. 
Similarly,  the  heat  dissipated  per  dollar  cost  of  the  earth  heat  sink 
for  the  2-week  and  3-week  duration  eart  heat  sinks  is  7200  Btu  and 
7700  Btu,  respectively.  Table  10  details  results  of  calculations  for 
the  two  types  of  soil  considered  here.  With  optimization  involving 
pipe  size,  pipe  spacing  and  the  material  for  select  backfill,  it 
should  not  be  difficult  to  increase  heat  dissipation  per  dollar  cost 
above  those  shown  in  Table  10.  Using  select  backfill,  the  heat 
dissipation  per  dollar  cost  of  between  5300  and  7700  Btu  appears  to 
be  relatively  inexpensive. 

Looking  at  it  from  another  viewpoint,  calculations  were  carried 
out  to  determine  the  Btu  absorption  capacity  of  the  earth  heat  sink  in 
Btu/cu  ft,  as  shown  in  Table  11.  If  water  at  60°F  or  ice  at  32°F 

j£o^lliHn4eKfirflurt.was.H^-  absor^tioi^capacities  will _ 

be  approximately  7,200  Btu  and  16,400  Btu  per  cubic  foot,  respecti^eryy 
Actual  volumetric  advantages  of  using  water  or  ice  will  be  somewhat 
less  than  that  shown  by  a  comparison  of  these  figures  with  those  of 
Table  11  since  the  volume  of  tank  to  hold  water  or  ice,  volume  of 
insulation  surrounding  the  ice  tank,  the  volume  of  refrigeration 
equipments  needed  to  produce  and  maintain  ice  and  the  volume  of  heat 
exchangers  needed  to  affect  the  heat  transfer  were  not  included  in 
these  estimates. 

Meaningful  cost  comparisons  with  alternative  methods  of  heat 
dissipation  require  a  knowledge  of  overpressures,  duration  of  use, 
cost  of  protecting  the  earth  heat  sink  against  ground  motion,  and  the 
cost  of  alternative  methods  of  heat  dissipation  and  their  development 
cost.  Most  of  this  information  is  not  available  at  this  time. 
Comparison  of  earth  heat  sinks  with  other  heat  dissipation  methods  are, 
however,  now  being  carried  out  for  conditions  applicable  to  project 
SANGUINE.  Preliminary  comparisons  are  encouraging. 

CONCLUSIONS 

Preliminary  calculations  have  shown  that  earth  heat  sinks  may 
be  practical  for  dissipation  of  waste  heat  from  underground  power 
plains.  Preliminary  calculations  show  that  an  earth  heat  sink  using 
select  backfill  may  dissipate  approximately  6000  to  8000  Btu  per  $ 
cost  and  5000  to  6000  Btu  per  cubic  foot  volume  of  the  earth  heat 
sink.  Optimization  of  design  may  produce  further  improvement  in 
performance  of  the  earth  heat  sinks. 
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RECOMMENDATIONS 


1.  A  complete  solution  should  be  obtained  for  the  case  of  a  long, 
concentric  circular  cylinder  whose  ends  and  outer  wall  are 
perfectly  insulated  and  whose  inner  wall  is  raised  and  maintained 

at  a  fixed  temperature  above  the  uniform,  initial  cylinder  temperature. 
The  solution  should  be  plotted  in  graphical  form  for  values  of  the  soil 
properties,  the  inner  and  outer  diameters  and  the  operating  period 
varying  over  wide  ranges. 

2.  Using  the  results  of  the  above  analysis,  an  optimization  of 
earth  heat  sinks  should  be  carried  out  to  obtain  the  lowest  overall 
volume  and  the  lowest  overall  cost. 


i 

j 


Table  1.  Thermal  Conductivities  of  Two 
Soil  Samples  at  60°F,  [28]  . 


Sample 

Density 
lb/cu  ft 

Moisture 

1 

k 

Btu/hr-f t-°F 

82 

9.9 

0.6 

86 

4.7 

0.3 

91 

9.9 

0.7 

Sandy- type 

95.6 

4.7 

0.45 

soil 

100 

9.9 

1.0 

103 

11.9 

1.25 

105 

4.7 

0.6 

112 

11.9 

1.47 

48 

41.4 

0.33 

51 

57 

0.43 

Clay  and  Loom- 

52 

29.7 

0.28 

type  soil 

53 

41.3 

0.37 

57 

57 

0.60 

59 

29.7 

0.33 

60 

41.3 

0.42 
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Table  2.  Location  of  various  pipes  in  a 
24  x  24  hexagonal  prism  matrix. 


Location 
of  pipes 
(Fig.  2) 

No.  of 
pipes  at 
similar 
locations 

Number  of  pipes  at  a 

distance  of 

2r 

0 

2VTr 

0 

4r 

0 

1 

2 

2 

2 

2 

2 

2 

4 

2 

3 

3 

40 

4 

3 

4 

4 

2 

5 

3 

2 

5 

2 

6 

5 

3 

6 

40 

6 

5 

4 

7 

20 

3 

3 

3 

8 

20 

6 

4 

4 

9 

400 

6 

6 

6 

10 

2 

4 

3 

3 

11 

2 

3 

2 

2 

12 

2 

3 

3 

2 

13 

2 

6 

4 

3 

14 

20 

5 

4 

3 

15 

20 

6 

6 

5 
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Table  3.  Values  of  I(Xn)  for  pipes  in  a 

24  x  24  hexagonal  prism  matrix 
in  an  earth  heat  sink  of 
Of  =  0.0185  sq  ft/hr. 


Pipe  Spacing 
(in) 

2v2TF 

(ft) 

r 

(ft) 

X 

n 

ny 

25 

3.52 

0.0625 

0.0178 

3.7780 

2.083 

0.5918 

0.3970 

3.610 

1.0256 

0.1007 

4.167 

1.1837 

0.0585 

30 

4.98 

0.0625 

0.0126 

4.1390 

2.50 

0.5020 

0.5191 

4.33 

0.8701 

0.1680 

5.00 

1.0040 

0.1083 

40 

6.10 

0.0625 

0.0103 

4.2990 

3.333 

0.5464 

0.4548 

5.770 

0.9459 

0.1313 

6.667 

1.0929 

0.0810 

Table  4.  Heat  transfer  rates  using  I(Xn) 

values  from  Table  3  for  pipe 
locations  shown  in  Figure  2  in 
an  earth  heat  sink  of  9  =  0.0185 
sq  ft/hr,  and  k  =  0.5  Btu/hr-ft-°F 


Location 
of  pipe 
in  the 
matrix 


SI(X  ) 
n 


Q'  (Btu/hr-f t) 


1-week  2-week  3-week  1-week  2-week  3-week 


4.8904 
5.7429 
5.9021 
6.1821 
6.8590 
6. 8975 
5.4466 
6.7968 
7.1152 
5.8436 
5.2874 
5.3881 
6.7383 
6.3413 
7.0567 


5.7298 

6.8763 

7.1526 

7.4551 

8.4135 

8.5268 

6.5282 

8.3588 

8.9114 

7.0443 

6.2489 

6.4169 

8.2505 

7.7314 

8.8031 


6.5332 

6.6238 

6.8361 

7.1289 

7.9273 

8.0083 

6.3003 

7.8770 

8.3016 

6.7551 

6.0880 

6.2193 

7.7960 

7.3412 

8.2206 


Weighted  average  for  the  24  x  24 
matrix - 


see  spacings  between  pipes  in  Table  3. 


Table  5. 


Dimensions  of  earth  heat  sinks  of 
a  =  0.0185  sq  ft/hr  and  k  =  0.5  Btu/hr-f t-°F 
for  a  300  kw  Diesel  using  a  24  x  24  hexagonal 
prism  pipe  matrix. 


Duration  of 
Operation 
(hr) 

168 


Pipe 

Spacing 

(in) 


Dimensions  of  the  Heat  Sink 


Length 

(ft) 

Width 

(ft) 

Depth 

(ft) 

Volume 
(cu  ft) 

91 

51 

44 

204000 

114 

61 

52 

362000 

106 

82 

70 

608000 

-  “  r’T  ’ 


Table  6.  Heat  transfer  rates  for  a  pipe  in  location  9 
in  the  Figure  2  matrix  in  an  earth  heat  sink 
of  a  =  0.0267  sq  ft/hr  and  k  =  1.2  Btu/hr-f t-°F. 


Time 

(hr) 

■ 

2I<V 

for 

location 

9 

Q' 

(Btu/hr-f t) 

168 

25 

4.2332 

0.0625 

0.0148 

3.9867 

8.983 

127.6 

2.083 

0.4921 

0.5395 

3.610 

0.8528 

0.1776 

4.167 

0.9844 

0.1156 

336 

30 

5.9867 

0.0625 

0.0104 

4.2861 

11.021 

104.0 

2.500 

0.4176 

0.6683 

4.333 

0.7233 

0.2663 

5.000 

0.8352 

0.1878 

504 

40 

7.3320 

0.0625 

0.0085 

4.4785 

10.267 

111.6 

3.333 

0.4546 

0.5982 

5.765 

0.7870 

0.2186 

6.667 

0.9093 

0.1480 

Table  7.  Dimensions  of  earth  heat  sinks  of  a  =  0.0267 
sq  ft/hr  and  k  =  1.2  Btu/hr-f t-°F  for  a  300 
kw  Diesel  using  a  24  x  24  pipe  matrix. 


Duration 

(hr) 

Pipe 

dia. 

(in) 

Pipe 

Spacing 

(in) 

Dimensions  of  the  heat  sink 

Length 

(ft) 

Width 

(ft) 

Depth 

(ft) 

Volume 
cu  ft 

168 

1.5 

25 

49 

51 

44 

110000 

336 

1.5 

30 

60 

61 

52 

190000 

504 

1.5 

40 

57 

82 

70 

327000 

26 


Table  8.  Heat  transfer  rates  for  a  pipe  in 

location  9  in  the  Figure  2  matrix  in 
an  earth  heat  sink  of  a  =  0.0267  sq  ft/hr 
and  k  =  1.2  Btu/hr-f t-°F. 


Time 

(hr) 

Pipe 

Spac¬ 

ing 

(in) 

2vtfF 

(ft) 

r 

(ft) 

X 

n 

2KX  ) 
for0 

location 

9 

Q' 

(Btu/hr-f t) 

168 

30 

4.2332 

0.0625 

0.0148 

3.9867 

7.3437 

156.0 

2.500 

0.5906 

0.3984 

4.325 

1.0217 

0.1020 

5.000 

1.1811 

0.0591 

336 

40 

5.9867 

0.0625 

0.0104 

4.2861 

8.1249 

141.0 

3.333 

0.5567 

0.4409 

5.765 

0.9630 

0.1241 

6.667 

1.1136 

0.0748 

504 

50 

7.3320 

0.0625 

0.0085 

4.4785 

8.1469 

140.7 

4.167 

0.5683 

0.4260 

7.210 

0.9834 

0.1160 

8.333 

1.1366 

0.0694 

Table  9.  Dimensions  of  earth  heat  sinks  of 
0t  =  0.0267  sq  ft/hr,  and  k  =  1.2 
Btu/hr-ft-°F  for  a  300  kw  Diesel  using 
a  24  x  24  pipe  matrix. 


Duration 

(hr) 


Pipe 

Pipe 

dia. 

Spacing 

(in) 

(in) 

Dimensions  of  the  heat  sink 


Length 

(ft) 

Width 

(ft) 

Depth 

(ft) 

41 

61 

50 

46 

82 

70 

47 

102 

87 

Vo lume 
cu  ft 


168 

336 

504 


30 

40 

50 


130000 

264000 

417000 
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NOMENCLATURE 


0^  Specific  heat  of  soil 


Btu/lb-°F 


Jq,  Bessel  functions  of  the  first  kind 

k  Thermal  Conductivity  of  soil 

n  Number  of  pipes  in  the  array 

N^  Fourier's  number,  r^/gt 

p  The  ratio  of  distance  of  a  location 
under  consideration  from  the  center 
of  the  pipe  to  the  radius  of  the  pipe. 

Q '  Heat  transfer  rate  per  foot  of  pipe 

r  Distance  from  center  of  pipe 

Radius  of  the  pipe 

r2>  r3»  etc  Distance  of  heat  sources  or  sinks  from 
the  center  of  pipe  of  radius  r^ 

r^,  Inner  radius  of  a  cylinder 

r  Outer  radius  of  a  cylinder 


Btu/hr-f t-°F 


Btu/hr-ft 


T  Temperature  at  distance  r  at  time  t 

T^  Temperature  of  pipe  of  radius  r^ 
at  time  t 

T  Equilibrium  earth  temperature  at 
time  zero 

t  Time 

X  r/2  \ST 


Y  ,  Y,  Bessel  function  of  the  second  kind 
o'  1 


A  W- 


